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Abstract 

Quantum Hamiltonians that are fine-tuned to their so-called Rokhsar-Kivelson (RK) points, 
first presented in the context of quantum dimer models, are defined by their representations 
in preferred bases in which their ground state wave functions are intimately related to the 
partition functions of combinatorial problems of classical statistical physics. We show that 
all the known examples of quantum Hamiltonians, when fine-tuned to their RK points, be- 
long to a larger class of real, symmetric, and irreducible matrices that admit what we dub 
a Stochastic Matrix Form (SMF) decomposition. Matrices that are SMF decomposable are 
shown to be in one-to-one correspondence with stochastic classical systems described by a 
Master equation of the matrix type, hence their name. It then follows that the equilibrium 
partition function of the stochastic classical system partly controls the zero-temperature 
quantum phase diagram, while the relaxation rates of the stochastic classical system co- 
incide with the excitation spectrum of the quantum problem. Given a generic quantum 
Hamiltonian construed as an abstract operator defined on some Hilbert space, we prove 
that there exists a continuous manifold of bases in which the representation of the quan- 
tum Hamiltonian is SMF decomposable, i.e., there is a (continuous) manifold of distinct 
stochastic classical systems related to the same quantum problem. Finally, we illustrate 
with three examples of Hamiltonians fine-tuned to their RK points, the triangular quan- 
tum dimer model, the quantum eight-vertex model, and the quantum three-coloring model 
on the honeycomb lattice, how they can be understood within our framework, and how 
this allows for immediate generalizations, e.g., by adding non-trivial interactions to these 
models. 
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1 Introduction 



The combinatorial problem of counting how many ways there are to pack dimers on 
somegiven lattice is of relevance to chemists, physicists, and mathematicians fH 
A deep connection between the statistical physics of closely-packed, hard- 
core dimers on two-dimensional lattices, and the critical behavior of the two-dimensional 
Ising model was established in the early 60's ||2L|3j]. Soon after the discovery of 
high-temperature superconductivity, a quantum version of the classical hard-core 
dimer problem on the square lattice was proposed by Kivelson, Rokhsar, and Sethna 



as an effective low energy theory for a doped Mott insulator M5Ll6U7ll8Ll9 Ul0l.il 111 . 



This example of a square lattice quantum dimer model is interesting and unusual 
in several ways. (1) There is a one-to-one correspondence between the (dimer) 
basis that spans the underlying Hilbert space and the configuration space of the 
combinatorial problem. (2) The quantum Hamiltonian is the sum over local Hermi- 
tian operators, each of which encodes the competition between a potential energy 
that favors a local ordering of dimers and a kinetic energy term that favors a lo- 
cal quantum superposition of dimers, i.e., a local quantum liquid state. (3) For a 
special value of the ratio between the characteristic potential and kinetic energies 
called the Rokhsar-Kivelson (RK) point, the local Hermitian operators entering the 
Hamiltonian are positive semidefinite and the ground state (GS) is the equal-weight 
superposition of all dimer states, i.e., the normalization of the GS is nothing but the 
number of ways to closely -pack hard-core dimers on the square lattice [6]. Remark- 
ably, this GS is a very peculiar liquid state since it is critical according to the results 
of Kasteleyn on the classical square lattice dimer model [2]. (4) It was realized by 
Henley that the excitation spectrum of the quantum square lattice dimer model at 
the RK point is identical to the spectrum of relaxation rates of the classical square 
lattice dimer model out of thermal equilibrium when equipped with a properly cho- 
sen Monte Carlo dynamics ill 211 . 

Building on the close interplay between the classical and quantum square lattice 
dimer models, Moessner and Sondhi extended the triangular lattice classical dimer 
model to a quantum one and showed that its RK point realizes an incompressible 
spin liquid state characterized by topological quantum order and robust to small 
perturbations iflill . The triangular lattice quantum dimer model provides a theoreti- 
cal playground in which concepts related to fractional quantum numbers, emergent 
gauge symmetries, and deconfinement of gauge quantum numbers take a precise 
form IU4I1 . This work triggered the study of classical and quantum dimer models 
on a variety of lattices in two and three dimensions, all quantum ones being char- 
acterized by RK points, but differing on t he p resence or not of a two-sublattice 
structure 0B|3HBBHSQ[21[23. 
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It was noted by Ardonne, Fendley, and Fradkin that quantum Hamiltonians that can 
be fine-tuned to an RK point are not exclusively built from classical dimer models 
and that the RK point can be extended to a line, or, more generally, to a higher- 
dimensional region of parameter space W26I1 . More specifically, quantum six- and 
eight- vertex models as well as quantum Hamiltonians not anymore related to com- 
binatorial statistical physics were shown to possess a manifold of RK points lE^l . 
Around the same time, Henley pointed out that any classical system endowed with 
linear dynamics in time (through a Master equation of the matrix type) can be used 
to construct a quantum Hamiltonian at an RK point H27I1 . 



In spite of the vast literature on the properties of quantum Hamiltonians fine-tuned 
to their RK points, there are still some fundamental questions that remain unan- 
swered. What are the generic properties of a quantum Hamiltonian fine-tuned to an 
RK point, and can they be used to define a more general class of problems? Under 
what conditions can one fine-tune the coupling constants of a quantum Hamiltonian 
to reach an RK point? Are all quantum Hamiltonians in this general class related 
to classical statistical systems and if so what is the precise nature (static versus 
dynamic) of this relationship? Given that the very notion of an RK point is basis- 
dependent, what is the fate of an RK point under a change of basis? In this paper 
we will address and answer these questions, thus providing insights into the re- 
markable relationship between quantum systems fine-tuned to their RK points and 
stochastic classical systems described by a Master equation of the matrix type. 

Common to all known quantum Hamiltonians when fine-tuned to their RK points 
are the following properties: 

(I) The underlying Hilbert H space is separable, i.e., the set S = {C} of indices 

C labeling a basis B = {\C), C G S] of TC is countable. 
(II) The quantum Hamiltonian can be decomposed into the sum of positive- 
semidefinite Hermitian operators, each of which being equal to its square after 
proper rescaling, weighted by positive coefficients. In particular, when writ- 
ten in the preferred basis, these operators are represented by 2 x 2 matrices 
acting solely on a two-dimensional subspace of the Hilbert space H. Any of 
these operators encodes the competition between a potential energy diagonal 
in the preferred basis and a kinetic energy term that favors a pairwise quantum 
superposition of basis states. 
(Ill) A GS of the quantum Hamiltonian can be found by demanding that it is 
annihilated by each and every Hermitian operator defined in (II). As we shall 
see, the normalization factor in the GS wavefunction can be interpreted as the 
partition function of a specific classical system defined over the phase space 
S. 



We shall first construct the most general representation of a quantum Hamiltonian, 
given a preferred basis B = {\C), C G S}, satisfying conditions (I-III), and we 
will show that it has a parametrical form that we shall dub the Stochastic Matrix 
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Form (SMF) decomposition. The set of all matrices that are SMF decomposable in- 
clude all known examples of quantum Hamiltonians fine-tuned to their RK points 
but are not limited to these. We shall show that conditions (II) and (III) are repre- 
sented by the positive-semidefinite decomposition conditions [Eq. (8)] and by the 
integrability conditions [Eq. (14)], respectively. In the preferred basis, the positive- 
semidefinite decomposition conditions have a local character, while the integra- 
bility conditions have a global character. Both set of conditions are automatically 
satisfied when quantum Hamiltonians fine-tuned to an RK point are constructed 
starting from a classical system as in Refs. [26] and lEvll . say. Whereas the im- 
portance of the local positive-semidefinite decomposition conditions has long been 
implicitly known, the fact that they do not necessarily imply the global integrabil- 
ity conditions has been overlooked so far. The global integrability conditions play 
a crucial role when the SMF decomposition of a quantum Hamiltonian represented 
in the preferred basis is to be achieved by fine-tuning of its coupling constants. 

In answering the third question, whether all representations of quantum Hamilto- 
nians that are SMF decomposable are related to classical statistical systems and if 
so what is the precise nature of this relationship, we will first show that it is al- 
ways possible to extract from the GS wavefunction of a quantum Hamiltonian that 
is SMF decomposable, when written in the preferred basis, the partition function 
of a classical system in thermal equilibrium. It then follows that the equal-time 
correlation functions (of operators diagonal in the preferred basis) in this GS can 
be interpreted as equal-time correlation functions in a classical system in thermal 
equilibrium. In that sense, the quantum phase diagram of a quantum Hamiltonian 
that is SMF decomposable when written in the preferred basis contains the thermal 
phase diagram of a classical system in thermal equilibrium. This classical system 
needs not be a pure combinatorial problem, nor one equipped with some chemical 
potential. Instead, it is generically described by a partition function evaluated at the 
inverse temperature (3 and given by a sum over Boltzmann weights with config- 
uration energies E c for all C E S. Moreover, as a direct consequence of proper- 
ties (II- III), we will show that the matrix elements of the quantum Hamiltonian in 
the preferred basis B = {\C), C E S} can be naturally interpreted, via a similarity 
transformation, as transition rates of a Master equation of the matrix type for the 
stochastic evolution of the corresponding classical system with configuration space 
S. These transition rates satisfy detailed balance precisely with respect to the equi- 
librium probability distribution induced by the GS wavefunction. Consequently, the 
excitation spectrum of any quantum problem that admits an SMF representation in 
the preferred basis coincides with the relaxation rates of a stochastic classical sys- 
tem and, as was observed by Henley [27], quantum correlation functions at unequal 
imaginary times in the GS for operators diagonal in the basis B — {\C), C E S} 
are identical to classical correlation functions at unequal times. We will show that 
quantum Hamiltonians that admit an SMF representation in the preferred basis are 
parametrized by the (independent) positive coupling constants w c c , and by the (in- 
dependent) coupling constants E c . Whereas the coupling constants E c determine 
the GS wavefunction, the non-vanishing coupling constants w c c , solely determine 
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Table 1 

Correspondences between the SMF representation of the quantum Hamiltonian H SMF 
from Eq. (15) and the classical system with stochastic dynamics in time captured by the 
transition matrix (W cc ,) from Eq. (19). 



Quantum system admitting an SMF decomposition Classical system 



Hilbert space with basis B labeled by S 
Ground state wavefunction 
Quantum phase transitions 
Hamiltonian matrix: (^(C\H SMF \C')^j 

Positive-semidefinite 

decomposition conditions 

Integrability conditions 

Energy eigenvalues 
Eigenfunctions 



Configuration space S 
Boltzmann distribution 
Classical phase transitions 
Transition matrix: (W cc ,) 

Positive transition rates 

Conservation of probabilities 

Relaxation rates 
Right/Left eigenfunctions 



the excitation spectrum above the GS wavefunction. Remarkably, the excitation 
spectrum can be manipulated by fine-tuning the w c c , without affecting the nature 
of the GS wavefunction! This corresponds, in the associated classical system, to 
being able to change the stochastic dynamics without affecting the equilibrium par- 
tition function. We will also establish that representations of quantum Hamiltoni- 
ans, in some preferred basis B — {\C), C E S}, that are SMF decomposable are in 
one-to-one correspondence with discrete classical systems with configuration space 
S endowed with stochastic dynamics described by a Master equation of the matrix 
type. The dictionary between the quantum and classical systems is summarized in 
Table 1. 



So far, our generalization of the known quantum Hamiltonians that are fine-tuned 
to their RK points is predicated on choosing a preferred basis first, i.e., given a 
special choice of basis, we have investigated under what conditions a quantum 
Hamiltonian is SMF decomposable. This is so for historical reasons. One can of 
course ask the reverse question, namely is there a basis in which some given quan- 
tum Hamiltonian is SMF decomposable? We shall give an affirmative answer to 
this question for any quantum Hamiltonian defined on a finite dimensional Hilbert 
space as long as it has not too many degenerate eigenvalues, up to some trivial shift 
of the energy spectrum. In fact, there are continuously many bases in which such 
a quantum Hamiltonian admits distinct representations that are SMF decompos- 
able. The correspondence between a quantum Hamiltonian, when understood as an 
abstract operator acting on some Hilbert space, and stochastic classical systems is 
thus one-to-(continuously) many. 
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As an application of these results, we shall illustrate a systematic procedure to 
include dimer-dimer, vertex-vertex, or color-color interactions in known quantum 
models that are tuned to their RK points, such as the quantum dimer, vertex, and 
color models, respectively. In particular, we discuss the procedure in greater de- 
tail for the case of the quantum three-coloring model. These extra interactions can 
lead to additional quantum phases and, for quantum Hamiltonians that are SMF 
decomposable, such quantum phases and quantum phase transitions can be partly 
understood in terms of purely classical phase transitions at thermal equilibrium. A 
second facet of these results is that the quantum Hamiltonians have simple repre- 
sentations as sums over "elementary" operators with a rather intuitive interpretation 
in terms of "elementary" moves in the associated classical configuration spaces. For 
example, in the quantum three-coloring model, we shall first identify "decorated" 
loops in the classical version of the model, from which we construct the quantum 
Hamiltonian as a sum over "elementary" operators associated to these "decorated" 
objects. Each of these "elementary" operators corresponds to an "elementary" move 
in the associated classical phase space that is energy non-conserving. A similar pro- 
cedure can be carried out in quantum dimer and vertex models after having identi- 
fied "decorated" plaquettes that encode non-trivial dimer-dimer and vertex-vertex 
interactions, respectively. 

The first main result of our paper, in the form of the generic conditions on a given 
representation of a quantum Hamiltonian for it to be SMF decomposable, is derived 
in Sec. 2. In Sec. 3 we show how the quantum dynamics of a quantum Hamiltonian 
that admits a representation that is SMF decomposable induce classical stochas- 
tic dynamics encoded by an appropriate Master equation. This result, when com- 
bined with Henley's result in the opposite direction, H27I1 allows us to establish the 
one-to-one correspondence between all possible SMF decompositions of quantum 
Hamiltonians in the preferred basis B = {\C), C G S} and discrete, stochastic 
classical systems described by a Master equation of the matrix type in the classical 
configuration space S. In Sec. 4 we present the second main result of the paper, 
namely the conditions under which a quantum Hamiltonian admits a continuous 
manifold of bases for which it is SMF decomposable. Specific examples of how 
known quantum Hamiltonians that are fine-tuned to their RK points can be under- 
stood within our framework and how this allows for immediate generalizations are 
given in Sec. 5. These examples are the triangular dimer model, the eight- vertex 
model, and the three-coloring model on the honeycomb lattice. Conclusions are 
drawn in Sec. 6. 



2 The Stochastic Matrix Form decomposition of quantum Hamiltonians 

In this section we are going to construct the most general form of a quantum Hamil- 
tonian, when represented in a specific basis, satisfying the properties (I-III). 
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Consider the Hilbert space H given by the span of all orthogonal and normalized 
(orthonormal) states \C), labeled by the index C E S with S a countable set (condi- 
tion (I)), 

(C\C') = 6 CC „ t = J2\C)(C\. (1) 

ces 

Consider a generic (Hermitian) Hamiltonian H acting on this Hilbert space, and 
define the subset C C S x S to be the set of all the pairs (C, C) with C ^ C such 
that the off-diagonal matrix elements of the |<S| x \S\ Hermitian matrix 

(H cc ,) := ((C\H\C)) (2) 

are non-vanishing. For simplicity, we will make the two technical assumptions that 
the Hilbert space is finite dimensional and fully connected under the time-evolution 
operator, i.e., any two states in TC have a non-vanishing matrix element for some 
power of H. These two conditions will be needed when using the Perron-Frobenius 
theorem to establish the non-degeneracy of the GS. Then, we can always represent 
# in the basis £ = {\C), C E S} by 0] 

H=\ E he ( 3a ) 

z (Cfi')ec 

where 

h c ,c> ■= a C ,c>\C)(C\ + <Xe',c\C'){C'\ - lc,M(C'\ - lc , fi \C'){C\ (3b) 

and 

a c c , G R, 7 CjC , = 7c>,c e c ( 3c ) 
(7c c , is the complex conjugate to 7 C c ,). Observe that the representation (3) is re- 
dundant. Indeed, the matrix elements of H in this preferred basis, 

{C,c')ec 

H CC = E "co VCG5, (4) 

Ces 

H C c = -7c,c, V(C,C')e£, (5) 
are unchanged under the transformations 

a cc -* ®c,c> = a c,c + fc,ch (6) 

provided 

(c,c')ec 

0= E fee, VCG5, (7) 
C'eS 

where the two summations, in Eqs. (4) and (7), are over the single index C E 
S subject to the constraint that (C,C) E C. This is so because our Hermitian 
matrix (H cc ,) has |«S| independent real- valued diagonal entries and \C\/2 non- 
vanishing independent off-diagonal complex- valued matrix elements. On the other 
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hand, there are \C\ independent a's and \C\/2 independent 7's. The representa- 
tion (3) is thus redundant, at least for \S\ > 2, since \C\ > 2(\S\ — 1) (see [ 28]). 
It will nevertheless be useful to identify the set of data (C, cx cc ,, 7 C c ,) with the 
Hamiltonian H. 

Notice the competition between the potential energy, the first two terms of Eq. (3b), 
and the kinetic energy, the last two terms of Eq. (3b), in each Hermitian operator 
h c c ,. In order to satisfy condition (II), we have to require that 

a c,c x a c,c = hc,c\ 2 and a c,c > (8) 

for all (C,C) E C, in which case we are left with \C\/2 > (\S\ — 1) independent 
a's. Indeed, one verifies that Eqs. (3b) and (8) imply 

h 2 c ,c = { a c,c + a c,c) h C fi>- ( 9 ) 

The first condition in Eq. (8) amounts to nothing else but imposing that the 2x2 
Hermitian matrix 

/ 



(10) 



a c,c 7c,c 
\lc,c a c,c, 

defined on the two-dimensional vector space spanned by \C) and \C), has vanishing 
determinant. This is a necessary and sufficient condition for a Hermitian 2x2 
matrix to be proportional to a projection operator. Together, the two equations in (8) 
form the positive-semidefinite decomposition conditions, that are necessary and 
sufficient for the Hamiltonian H to be decomposable in 2 x 2 Hermitian blocks, each 
of which has precisely one null and one positive eigenvalue. Equation (8) is also 
sufficient to guarantee that H is a positive semidefinite quantum Hamiltonian. Any 
eigenstate annihilated by H thus necessarily belongs to the GS manifold whenever 
conditions (8) hold. 

Before imposing condition (III), it is useful to change the parametrization in Eq. (3b) 
to 



w c.c> '■= \ a c.c> x «c.c> ( lla ) 



e- K£ c^ := (Hb) 
V a c,c 

K<f>c,c/ 2 ■= arg(Tc,c)- ( llc ) 
Here, the coupling constants Ke c c , and K(p c c , are real valued and the w c c , are 
strictly positive. It is useful, even though redundant, to represent the real coupling 
constants Ke c c , and Kcf) c c , as products of a common scaling factor K and two 
state-dependent functions e c c , and C c ,. The reasons for this choice will become 
clear shortly, and rely on the physical interpretation of the scaling factor K (see 
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Eq. (15e)). After this change in parametrization, Eq. (3) becomes 

H=- w c,cQc,c ( 12a ) 

/ (c,c)ec 

where 

Q cc , = e- K£ cc' /2 \C)(C\ + e +K£ c-c l2 \C')(C'\ 
- e~ lK(i, e,ci/ 2 \C)(C\ - e +lK ^c>/ 2 \C){C'\. 



(12b) 



We now require that condition (III) be satisfied, i.e., that there exists a simultaneous 
zero mode for all the Hermitian operators Q c c , in Eq. (12b). One verifies that the 
nodeless wavefunction 

\%}:=Y.e- KEc/2 \C)> d3) 
ces 

is annihilated by Q cc ,, V (C, C) £ C, provided the integrability conditions 

£c,c = E c> ~ E c ( 14a ) 

and 

4> cc , = (up to a pure gauge transformation) (14b) 



on the real-valued parameters e c c , and <p c c , are satisfied for some real-valued func- 
tion E c defined on S, and for all (C, C) E C As a result of the integrability con- 
ditions (14b) and noting that all the couplings w cc , are positive by definition, we 

see that all off-diagonal matrix elements of H are negative. This property, together 
with the integrability conditions (14a) and the positive-semidefinite decomposi- 
tion conditions (8), will be needed to infer that H admits an SMF decomposition. 
Here, it suffices to say that the two conditions (14) guarantee that |\& ) is the GS 
of the Hamiltonian H in Eq. (12a), regardless of the values taken by the coupling 
constants w c c ,, and that this GS is non-degenerate by the Perron-Frobenius theo- 
rem H. 

We conclude that the most general form for a quantum Hamiltonian that is defined 
on the Hilbert space spanned by the basis B = {\C), C £ S} and that satisfies 
the RK conditions (I- III) is to depend on the data (C,w cc ,,E c ,K) through the 
representation 

Hsmf = 9 w c,cQc,c (15a) 

(c,c)ec 

where 

Q cc , = e - K{E c>- E cy 2 \C)(C\ + e +K{E c>- E c)l 2 \C')(C'\ - \C')(C\ - \C){C'\. (15b) 
This representation will be referred to hereafter as the Stochastic Matrix Form 
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(SMF) decomposition. The normalized GS wavefunction 



\^smf) = —n=7 E e~ KEc/2 \C) (15c) 
'Z(K) cg5 



of H SMF is non-degenerate as long as the Hilbert space is connected under the 
time-evolution operator and finite-dimensional. Here, the normalization factor 

Z(K) := ]T e~ KE c (15d) 
ees 

can be interpreted as the partition function of a classical system with phase space 
S, reduced inverse temperature 

K := (3 J, (15e) 

(J being some characteristic energy scale), and dimensionless energies E c . Evi- 
dently, in the infinite temperature limit (f3 = 0) and in a finite system, Z(0) reduces 
to the combinatorial problem of counting the number of classical configurations 
belonging to S. 

It is instructive to observe that Eq. (4) now takes the form 

(c,c')ec 



H cc = E 



w c c , e 



~ Ks c,c'/ 2 



C'&S 
(c,C)ec 

= E \H cc ,\e' Ke ^l\ (16a) 
ces 



(c,C)ec 



= E \Hcc>\^> VCG5, 
where we have introduced the notation 

e -KE c /2 

:= (C\* ) = —j==. (16b) 

As we shall prove in Sec. 4, for any quantum Hamiltonian H that has a ground 
state |\& ) with vanishing eigenvalue, and that admits an irreducible representation 
(H cc ,) with non-positive off-diagonal matrix elements in some finite dimensional 
basis B = {\C), C £ S}, it is possible to define in a unique way the decom- 
position (12) with all cf) cc , = 0, by defining the sets {Ke cc ,, (C,C) G £} and 
{ w cci e ^} through Ke cc ,/ 2 := ln(^ oc ,/^ oc ) and through w cc , : = 

\H CC ,\, so that Eq. (16a) is satisfied. We are now equipped with necessary and suf- 
ficient conditions for the matrix representation of a quantum Hamiltonian to be 
uniquely decomposable in the SMF (15). 
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Equations (15) encode our first important result: Given any quantum problem and 
a preferred basis, i.e., given a Hilbert space Ti. and a set of data (£, a c c ,, 7 C c ,) 
as defined above, it satisfies the RK conditions (I-III) iff Ti is separable and its 
coupling constants a c c , and 7 C c , satisfy the positive-semidefinite decomposition 
conditions (8) and the integrability conditions (14). The coupling constants can 
then be redefined 

(A a C,C'i 7c,c) (A w C,C'i E c K ) ( 17 ) 
and the Hamiltonian representation takes the simpler form (15). One can associate 
this quantum problem to a classical problem in thermal equilibrium described by 
the partition function (15d) appearing in the GS wavefunction (15c) of H SMF . All 
equal-time quantum correlation functions for operators diagonal in the preferred 
basis (1) of H are identical to their classical counterparts in thermal equilibrium at 
reduced inverse temperature K. Any finite temperature phase transition, say at K c , 
of the classical model (15d) implies a quantum (zero-temperature) phase transition 
at the RK point parametrized by the quantum coupling constant K c (the converse 
statement is not true!). Another remarkable property of H SMF is that, while the GS 
wavefunction is completely independent of the choice of the parameters w c c ,, they 
instead control to a large extent the nature of the excitations above the GS (15c). 

Conversely, for any discrete classical statistical system in thermal equilibrium as 
defined by the partition function (15d), say, we can construct the Hamiltonian 
H SMF with GS wavefunction (15c). Of course, H SMF is not uniquely defined as 
the coupling constants w c c , in the set of data (£, w c c ,, E c , K) needed to specify 

uniquely H SMF cannot be extracted from the equilibrium properties of the classical 
system. For the full correspondence between quantum Hamiltonian representations 
that are SMF decomposable and classical statistical systems to be established, we 
need also to account for the approach to thermal equilibrium in the classical system 
as we now explain. 



3 From quantum dynamics to classical stochastics 

Henley observed that the eigenvalues of the square lattice quantum dimer Hamil- 
tonian at its RK point are identical to the relaxation rates of a (classical) Mas- 
ter equation 11211 . Ivanov took advantage of this observation to simulate the vison 
gap of the triangular lattice quantum dimer Hamiltonian at its RK point through a 
classical Monte Carlo simulation on the same lattice as the one for the quantum 
model 13 ill . Henley also observed that any generic Master equation can be brought 
to the form of a quantum Hamiltonian fine-tuned to its RK point through a simi- 
larity transformation in a unique fashion [27]. The reverse question of whether any 
quantum Hamiltonian, when fine-tuned to an RK point, can be brought to the form 
of a Master equation has not yet been addressed. 

In this section we will construct a (classical) Master equation from any given 
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quantum Hamiltonian that admits an SMF decomposition when represented in a 
preferred basis, such that its stationary probability distribution is nothing but the 
Boltzmann weights from Eq. (15d). This Master equation can be thought of as the 
result of an environment that endows the classical system with specific stochas- 
tics in time. We will then show that the positive-semidefinite decomposition condi- 
tions (8) and the integrability conditions (14) combine in such a way that the Master 
equation can be written in matrix form, and the eigenvalues of the transition ma- 
trix (i.e., the relaxation rates of the classical system) are identical, up to an overall 
sign, to the positive eigenvalues of the quantum Hamiltonian (15). Here, the (local) 
positive-semidefinite decomposition conditions (8) are sufficient for all eigenvalues 
of the transition matrix to be non-positive, and together with the (global) integra- 
bility conditions (14) they are necessary for the transition rates to be non-negative 
and for the eigenstates of the transition matrix to have a conserved normalization 
under time evolution (i.e., conservation of probability holds). This establishes the 
correspondence between quantum dynamics and classical stochastics in the com- 
plementary direction to the one explored by Henley. The general character of the 
classical systems associated to SMF representations of quantum Hamiltonians al- 
lows us to present the following result: There exists a one-to-one correspondence 
between quantum Hamiltonians that admit an SMF decomposition (15) in a given 
preferred basis B = {\C), C E S} and classical statistical systems endowed with 
time stochastics through a Master equation of the matrix type (22) in the given 
configuration space S. 



3. 1 Master equation and transition matrix 



From now on, we denote the eigenvalues and the orthonormal eigenstates of H SMF 
in Eq. (15) by e n and \e n ) respectively, i.e., 

HsMF = ^2 £ n \ £ n) ( £ n\ (18) 
n 

with the non-degenerate GS |£ n=0 ) = \^smf)- Choose any two C and C in S and 
define the matrix elements 

H cc , := (C\H SMF \C) (19a) 

and 

W cc , := -e- K ( E c- E o)/ 2 H cc „ (19b) 

respectively. Hermiticity and time-reversal symmetry of H SMF imply the condition 
of symmetry 

H CC , = H C , C (20a) 
which, in turn, implies the condition of detailed balance 

W cc ,e- KE o = W c , c e~ KE c (20b) 
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for any pair C, C G S. 



Because of the integrability conditions (14b), given Eqs. (1 1,12) and definitions (19), 
all off-diagonal matrix elements in (W cc ,) are positive and can thus be interpreted 
as transition rates (conditional probabilities). We are now in the position to define 
the Master equation 



c'+c 

Pc( T ) ■= E W cc , Pc/ (t) - W c , cPc (t) 
Ces 



VCeS (21) 



whose properly normalized solution can be interpreted as the instantaneous prob- 
ability (a number between and 1) for the classical system to be in configuration 
C. Equation (21) defines in a natural and unique way the classical stochastics at 
the reduced temperature K induced by the quantum Hamiltonian (15) on the as- 
sociated classical system in thermal equilibrium. On the other hand, the special 
balance between the kinetic and the potential terms characteristic of a quantum 
Hamiltonian that admits an SMF decomposition when represented in a preferred 
basis guarantees that the transition matrix (W cc ,) satisfies the conservation of nor- 
malization condition W cc = — J2ces W c , c , VC G S, as it can be verified directly 
using Eqs. (15a,15b) and (19b). This is indeed a condition relating the diagonal 
and off-diagonal elements of (W cc ,) which follows from the SMF conditions (8) 
and (14). We are then allowed to recast Eq. (21) as a Master equation of the matrix 
type 

Pc(r) = E W cc , Pc ,(t), VC G S, (22a) 
Ces 

where the (transition matrix) positivity conditions 

W cc , > VC, C G S, C ^ C, (22b) 

the conservation of probability conditions 

eve 

W cc = ~ E W cc VC G S, (22c) 

ces 

and the detailed balance conditions 

Way pf = W c , c pf V C, C 6 5, (22d) 



^'-zk e ' KEc - <22e) 

hold. The solutions to Eq. (22) can be interpreted as probabilities at all times given 
that they are probabilities initially. They are of the form 

Pc(T) = E a ne- lKlT 4 R '' n \ (23) 
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and relax to the Boltzmann distribution (22e) as r — > oo. Here, the right-eigenvalues 
{A n } and right-eigenvectors {ipc of the transition matrix (W cc ,) control the 
relaxational time-dependence of the solutions (23) to the Master equation (22). 
Note that the spectral representation of the (not necessarily symmetric) matrix 
(Wcc) requires the introduction of left-eigenvectors {V4 } m addition to the 
right-eigenvectors and takes the form 

W cc , = Y.K4 R;n) 4''' n \ (24) 

nGN 

^„ = E# ;n Vf m) - (25) 
ees 

Comparing Eq. (24) with the definition (19b), one can verify the one-to-one corre- 
spondence 

K = ~ £ m (26a) 

4 R '' n) = e- KE ^x(C\e n ), (26b) 

4^ = e +KE c ,2 x ^ (2fic) 

with the eigenvectors and eigenvalues of the quantum Hamiltonian H SMF in Eq. (18), 
in the preferred basis. This fundamental property is a direct consequence of the 
SMF decomposition of the Hamiltonian (22a) and plays a key role in establishing 
the equivalence between classical and quantum correlation functions, proved by 
Henley in Ref. H- 

We close the discussion on the classical stochastics induced by a quantum Hamilto- 
nian that admits an SMF decomposition by giving the explicit form of the transition 
matrix in the preferred basis. Substituting Eqs. (15a,15b) into the definitions (19) 
we obtain 

'^,ce-^-^ )/2 ,if(C,C')G£, 



W, 



CO 



C"+C 

E W. 



C"Ci 



if (C, C) £ £, 



ifC = C 



(27) 



Notice that the last line of Eq. (27) is guaranteed to hold by conditions (22c). For 
the purpose of performing numerical simulations of the classical system, it is worth 
noticing that the special cases of Metropolis and Glauber dynamics can be imple- 
mented by choosing 



(28a) 
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and 

W C,C> OC e _ K{Ec _ Ec)/2 ~ e+K (E c -E c )/2> < 28b) 

respectively (recall that w c c , is required to be symmetric upon exchanging C and 
C by construction). 

We have shown in Sec. 2 that the coefficients w c c , remain free after relating the 
GS of a quantum Hamiltonian that admits an SMF decomposition to the partition 
function of a classical system in thermal equilibrium. We now see that the coeffi- 
cients w c c , play a crucial role when extending the relation between the quantum 
and classical systems to the excitation and relaxation spectra, respectively. We note 
that there are precisely as many independent parameters w c c , as there are indepen- 
dent non- vanishing matrix elements in the transition matrix (W c c ,). This fact has 
two consequences: 

(1) If we equip the classical system (15d) with some stochastics that can be written as 
the Master equation of the matrix type (22), then there exists a unique set {w c c ,} 
that implements the quantum dynamics in some representation of the quantum 
Hamiltonian (£, w c c ,, E c , K) in the sense that relaxation and quantum normal 
modes are identical. Thus, the representation (15) for a quantum Hamiltonian 
is indeed the most general one that follows from using Henley's procedure in 
Ref. H to construct a quantum Hamiltonian satisfying the RK conditions (I- 
III) starting from a Master equation of the matrix type (22). 

(2) Starting from a quantum Hamiltonian that admits an SMF decomposition (£, w c c , , E c 
we have constructed a unique and generic classical system equipped with time 
stochastics through the Master equation (22) whose relaxation rates are identi- 
cal, up to an overall sign, to the energy eigenvalues of the quantum Hamiltonian 
(£,w CiC ,,E c ,K). 

Combining (1) and (2), we have established a one-to-one correspondence between 
quantum Hamiltonians that admit SMF decompositions (15) and stochastic classi- 
cal systems that can be represented by the Master equation of the matrix type (22). 



4 Two Theorems about the Stochastic Matrix Form decomposition 

In the two previous sections we have shown how one can deduce from the known 
examples of quantum Hamiltonians fine-tuned to their RK points a general class 
of matrix representations of quantum Hamiltonians that exhibit the same charac- 
teristic properties, and how these matrix representations are related to stochastic 
classical systems. In this construction, the positive semi-definite decomposition 
conditions (8) and the integrability conditions (14) play a crucial role. Both con- 
ditions depend not only on the choice of the preferred basis B = {\C), C E S} 
but also on the non-unique decomposition (4) of the diagonal matrix elements of 
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the Hamiltonian H. We shall show in this section under what conditions there ex- 
ists a unique decomposition of the diagonal elements of a quantum Hamiltonian 
if in a suitable irreducible representation that satisfies the positive semi-definite 
decomposition conditions (8) and the integrability conditions (14). 

Also, we have shown in Sec. 3 that, to any SMF representation of a quantum Hamil- 
tonian H, in some preferred basis B — {\C}, C E S}, there corresponds a unique 
Master equation of the matrix type. We shall give below sufficient conditions under 
which a quantum Hamiltonian H, when construed as an abstract operator acting 
on a Hilbert space, admits different SMF representations and thus corresponds to 
different classical systems described by Master equations of the matrix type. 

Theorem 1 Given any real Hamiltonian H defined on a finite dimensional Hilbert 
space 7i such that ( i) it has a vanishing ground state energy, ( ii) there exists a basis 
B = {\C), C G S} such that H is represented by the irreducible matrix (H cc ,), 
and (Hi) H cc , < forC ^ C, then (H cc ,) is an SMF decomposable representation 
which is unique in every such basis B, up to a permutation of the basis elements. 

An immediate corollary is that the RK conditions (I-III) in Sec. 1 are always sat- 
isfied by a matrix that meets the hypotheses of Theorem 1 . In this sense, there is 
nothing special about the known examples of quantum Hamiltonians tuned to their 
RK points with respect to any other value of the parameters, other than the fact that 
the SMF decomposition is promptly handed-in from the knowledge of the energies 
{E c , C G S} (through the ground state wavefunction) and the off-diagonal matrix 
elements H cc , < for C ^ C. As an example, we can apply Theorem 1 to the 
quantum dimer problem on the square lattice 0,|6J]. We recall that the quantum 
dimer Hamiltonian is parametrized by two coupling constants commonly denoted 
v G (— oo, +oo) and t G (0, oo) for the diagonal and off-diagonal contributions 
to the quantum dimer Hamiltonian represented in the usual classical-configuration- 
state basis, respectively. The RK point t = v is the one for which the ground state 
wave function can be constructed explicitly and shown to have a vanishing en- 
ergy eigenvalue. However, since the quantum dimer Hamiltonian is real- valued and 
its off-diagonal elements are non-positive, with the restriction to limit the Hilbert 
space to one of the sectors in which the Hamiltonian is irreducible it is possible to 
apply Theorem 1 for any value of the dimensionless coupling t/v up to a shift of 
the spectrum that ensures that the ground state has a vanishing eigenvalue. Hence, 
the quantum dimer Hamiltonian allows for an SMF decomposition in the classical 
dimer configuration basis for any value of t/v. In other words, it is always possi- 
ble to write the quantum dimer Hamiltonian as a sum of 2 x 2 projectors in the 
classical-configuration-state basis, but the zero-energy GS becomes a superposi- 
tion of classical configuration states with positive coefficients that are no longer 
required to be all equal. If so, the spectral properties of the quantum dimer model 
at any t/v other than the RK point t = v can also be recast as a problem in clas- 
sical stochastic dynamics and thus amenable to classical simulations as opposed to 
quantum simulations. Of course, aside from the problem of solving for the ground 
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state, it is likely that the classical stochastic dynamics cannot be implemented by 
local updates of the classical configurations induced by the SMF representation. 

An immediate generalization of Theorem 1 can be achieved by considering com- 
plex Hamiltonians that satisfy hypotheses (i) and (ii) but fail to satisfy hypothe- 
sis (iii), provided all the off-diagonal elements can be made real and non-positive 
via a gauge transformation, i.e., provided there exists a unitary diagonal basis trans- 
formation 

|C)-e^|C), c G(O,27r) (29) 
that makes all transformed off-diagonal matrix elements non-positive. 

On completely general grounds, one can formulate an even stronger theorem. 

Theorem 2 Any Hamiltonian H defined on a finite dimensional Hilbert space 7i 
with a non- degenerate ground state whose eigenvalue is vanishing admits a con- 
tinuous manifold of distinct block diagonal representations where each irreducible 
block is an independent SMF decomposition. 

Theorem 2 immediately implies that the correspondence between a quantum Hamil- 
tonian H satisfying the hypotheses of Theorem 2 and stochastic classical systems 
is one-to-(continuously) many. The assumption that the ground state has a vanish- 
ing eigenvalue is not essential since one can always achieve this by shifting rigidly 
the energy spectrum by the subtraction of the unit matrix multiplied by the ground 
state energy. 

We now turn to the proofs of Theorems 1 and 2. 



Proof (Theorem 1) Let (H cc ,) be the matrix representation of the quantum Hamil- 
tonian H that satisfies the hypotheses of Theorem 1. All the off-diagonal matrix el- 
ements of — (H cc ,) are non-negative. By choosing a suitable positive constant c, all 
matrix elements of —(H' ca ) := —(H cc ,) + c (5 CC ,) can be made non-negative. As 
such, the matrix —(H' cc ,) obeys the hypotheses of the strong form of the Perron- 
Frobenius theorem from Ref. W2911 (if the matrix —(H' cc ,) were to be reducible, 
i.e., (H cc ,) is reducible, it would obey only the hypotheses of the weak form 
of the Perron-Frobenius theorem from Ref. [29]. However, one can always work 
separately within each irreducible block). The application of the strong form of 
the Perron-Frobenius theorem to the matrix —(H' cc ,) implies that its eigenvector 
\[/q with the largest eigenvalue is a (non-degenerate) strictly positive eigenvector, 
i.e., \P 0C > 0, VC Since —(H' cc ,) and (H cc ,) share the same eigenvectors, the 
strictly positive eigenvector \l/ is also the ground state with vanishing eigenvalue 
of(H cc ,). 

We can then identify a set of Boltzmann weights from the strictly positive coeffi- 
cients of \l/ , assumed normalized to one, and we can read out the classical energies 
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{KE C := —2 In (*o,c) > ^ e ^} f rom these Boltzmann weights, apart from an 
irrelevant shift of the energies {E c , C E S} by the constant In [Z(K)]/K. Using 
these classical energies {E c , C E S} we can perform a similarity transformation 
on the matrix —(H cc ,) as was done in Sec. 3 and get the matrix (W cc ,) such that 
(W cc ,) has the left-eigenvector (1, 1, 1) with vanishing left-eigenvalue. As the 
off-diagonal elements of the matrix —(H cc ,) are non-negative and the Boltzmann 
weights are positive (\l/ is positive), the off-diagonal elements of (W cc ,) are also 
non-negative. Moreover, the vanishing left-eigenvalue of (1, 1, 1) ensures that 
the diagonal elements of (W cc ,) are minus the sum of the off-diagonal elements in 
the column to which they belong. Put together, this means that elements of (W cc ,) 
have a transition probability interpretation. If so we can make use of the fact that to 
each Master equation of the matrix type there corresponds a unique matrix obeying 
the SMF decomposition to conclude that the matrix (H cc ,) realizes an RK point. 

Since the components of the wave function \l/ are unique up to a permutation of 
the elements in the set S, so are the Boltzmann weights, and so is the SMF (15) of 
the matrix (H cc ,). □ 



Proof (Theorem 2) Let Hbea quantum Hamiltonian that satisfies the hypotheses 
of Theorem 2. We are going to construct a continuous manifold of matrix represen- 
tations of H that satisfy the hypotheses of Theorem 1 . 

We choose a basis of H that diagonalizes H, i.e., H admits the matrix representa- 
tion 

//; ; "A,- i,j = 0,...,|5|-l, (30) 
such that all its eigenvalues are ordered according to 

e < e i < ■■■ < e \s\-i- (31) 

By hypothesis, e = 0. Let A be any operator defined on H such that it can be 
represented by a real- valued antisymmetric matrix whose matrix elements 

A l3 = -A 3l (32) 

are all infinitesimal and strictly positive if % < j, 

Kj^ A l3 > 0. (33) 

Define 

:= e"%e +2 , S^H := e^He^ - H. (34) 

The Hamiltonian shares the same eigenvalues as H but is not diagonal in the 
basis that diagonalizes H. Because A is infinitesimal we can write 

5 (A) H=[H,A}. (35) 
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The matrix elements of S^H in the basis that diagonalizes H are 



(^) y = (e* - = 



(36) 



In view of Eqs. (31), (32), and (33), 



(5^H) tJ < 



(37) 



where the equality holds if there are degenerate eigenvalues or if % = j. We con- 
clude that the matrix representation 



obeys the hypotheses of Theorem 1 at least within every irreducible block, once the 
basis vectors are premuted in such a way as to make the representation block 
diagonal. Note that if none of the eigenstates has degeneracy larger than or equal to 
\S\/2, irreducibility is ensured. □ 



5 Examples 

Before reviewing a few examples of quantum Hamiltonians that admit SMF rep- 
resentations, it is worth mentioning that the notation chosen in Eqs. (15), although 
the most general, can usually be further simplified depending on the specific details 
of the system. In fact, independently of whether the starting point is a classical or 
quantum system, the set of changes in phase space (C, C) produced by the clas- 
sical stochastics or by the kinetic energy operator can often (although not always) 
be interpreted as local (on some lattice) rearrangements £ of degrees of freedom 
(e.g., spin flips or exchanges, plaquette dimer operations, loop updates). Under the 
condition that the energy change is local, all the relevant physics in a pair (C, C) 
depends only on a neighborhood of the local rearrangement I that connects C to C 
(i.e., E c , — E c = e t and w ce , = w e ). Then, the summation over pairs (C, C) G C 
can be reduced to a summation over local rearrangements t E C r , and the quan- 
tum Hamiltonian at its RK point assumes a simpler, more intuitive form. Practical 
examples of this simplification are given hereafter. 

We will first provide two known examples of Hamiltonians fine-tuned to their RK 
points and show how they relate to our general formalism. We will also suggest pos- 
sible generalizations of these models that arise naturally from our results. Then, we 
will provide an example for possible practical applications of the results presented 
in Sections 2 and 3. 



(38) 
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5. 1 The triangular lattice quantum dimer model 



The quantum dimer model on the triangular lattice is defined on the separable 
Hilbert space spanned by the preferred, orthonormal basis B = {\C), C E S}, 
where C stands for any of the possible close-packed, hardcore dimer coverings of 
the triangular lattice. This Hilbert space is not connected under the time-evolution 
operator generated by 

H := -tf + vV 

n p 3 (39) 

= EEH (ImXv^l + h.c.) + v chattel + , 

where the sum on % runs over all of the iV p plaquettes y\, making up the triangular 
lattice while the sum on a runs over all three even permutations of the three sublat- 



tices making up the triangular lattice 11311 . We call a plaquette occupied by a pair 



of parallel dimers flippable. This Hamiltonian is known to exhibit an RK point at 
v — t (— 1) lllill . At this RK point, one of the possible GS is the equal- weight su- 
perposition of all the elements of the preferred basis 11311 . Hence, the normalization 
of this GS is nothing but the purely combinatorial, classical problem of counting 
all the arrangements of hardcore and close-packed dimers on a triangular lattice. 

We are now going to show how Eq. (39) with v — t — 1 relates to Eq. (15) with 
K = 0, w c c , = 1, and S defined as in the last paragraph. First, we need to construct 
the set £ C SxS. To this end, we define C to be the set of all pairs of configurations 
in S x S that are mapped one onto the other by a single flippable plaquette update, 
i.e., C differs from C by one and only one flippable plaquette, say ^ for C and 
for C, whenever (C, C) E C. Equipped with C, we can then consider Eq. (15) 
with K = and w c c , = 1, and rearrange it using the fact that the operators \C)(C\ 
and \C) (C\ both enter Eq. (15) with the same weights. For any given plaquette 
and up to an even permutation of the three sublattices of the triangular lattice, it is 
possible to perform the two sums 1E2I1 

ImXml := E \C)(C\, (40a) 



and 



c > c' 

hA>te|:= E (40b) 



Here, ^ E C is a shortcut for {C, s.t. E C}, while C — > C stands for 
{C, s.t. C transforms into C upon the local rearrangement — > ^}. We have 
thus recovered Eq. (39) starting from Eq. (15) with K = and w cc , = 1. Of 
course there is nothing sacred with these two conditions. We might as well relax 
them without spoiling the RK properties (ITU) of our quantum Hamiltonian which 
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now takes the more general form 13311 



H S MF =2 E ™m [- (I^X^I + H - C -) 

rz±ec r (41) 

+e- X£ -/ 2 |^)(TA| + e +Xe -/ 2 |V\)(NA| 

whereby we have replaced the double sum in Eq. (39) with the single sum labeled 
by G C r , where C is the set of all possible plaquettes with two parallel dimers 
(i.e., or \^ for all positions i and directions a). The reasons for such a change 
in notation will be clear from the examples presented in the following Sections. 
Notice that a factor of 1/2 is now required to avoid double counting due to the 
distinction between the two dimer coverings and \A of the same plaquette. The 
Hamiltonian H SMF differs from the quantum Hamiltonian (39) fine-tuned to its RK 
point t = v by the presence of the classical interaction between the dimers in the 
form of the configuration energy E c , whereby E c , — E c = [34]. One of the 
possible GS of Eq. (41) has the classical partition function 

Z(K):=J2 ^ KE ° (42) 
ces 

for normalization. Hence, the phase diagram of this classical partition function 
along the reduced temperature axis K partly controls the zero-temperature phase 
diagram of the quantum Hamiltonian (41) as a function of K G R. Finally, the 
energy spectrum can be tuned by the coupling constants w c c , = w^. 



5.2 The quantum eight-vertex model 



We turn our attention to the quantum eight-vertex model at its RK point, as dis- 
cussed by Ardonne, Fendley, and Fradkin lE^l . The classical configuration space S 
consists of all possible arrangements of vertices of the eight types shown in Fig. 1 
so as to form a square lattice. The Boltzmann weight of a configuration C G S is 
obtained by taking the product of the positive numbers (fugacities) a 2 , b 2 , c 2 , and 
d 2 associated to all the vertices appearing in C, according to the pairings implied by 
Fig. 1 . We demand that the partition function respects the symmetries of the square 




aabbccdd 



Fig. 1. The allowed eight type of vertices with the corresponding real- valued numbers a, 
b, c, and d. 

lattice, in which case we must impose the condition a = b. For simplicity we set 
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a = b = 1. This gives the classical partition function 

Z(c 2 ,d 2 ):=Y,c 2N c C) d 2N c d) (43) 

where is the total number of vertices of type c, say, in configuration C. The 
separable Hilbert space H is the span of the preferred, orthonormal basis B = 
C G S}. In this basis, the quantum eight- vertex Hamiltonian fine-tuned to its 
RK point is represented by 

H RK :=J2 w iQn ( 44a ) 



Qi ■■= 1 , (44b) 




_(c) (c) —(d) (d) 
„, . „n —n jn: — n. r A A „\ 

v i := c 1 1 a 1 1 . (44c) 

Here, the index i, which runs over all the sites of the square lattice, labels all the 
elementary plaquettes of the square lattice (square plaquettes) through their lower- 
left corners, the coupling constant w { is a positive weight, and the 2x2 matrix 
Q i oc Q 2 relates any configuration C to the (unique) configuration C obtained from 
C by reversing all the arrows on the edges of plaquette i 13011 . The positive integers 
nf' and nf' count the numbers of vertices of type c, say, that appear on plaquette i 
in configuration C and C, respectively. 

The Hamiltonian (44) has the normalized GS wavefunction 

^Z(c 2 ,d 2 ) ces 

which is non-degenerate in any of its four topological sectors if periodic boundary 
conditions are imposed owing to the fact that all configurations are flippable lEfill . 
The coupling space of the quantum eight- vertex model fine-tuned to its RK point is 
two-dimensional and is parametrized by c and d, which can be taken to be positive 
without loss of generality. In this parameter space the quantum eight-vertex model 
fine-tuned to its RK point exhibits quantum phase transitions between a confining 
ordered phase, with broken Z 2 symmetry, and a deconfining disordered phase, with 
unbroken Z 2 symmetry 112611 . The lines of critical points at the boundaries between 
those two phases correspond to the classical six- vertex model case and its dual lE^l . 

In order to recast Eq. (44) into the formalism of the past sections, observe that the 
potential energy term acts multiplicatively on a basis state \C) with the proportion- 
ality constant given by 

_(c) (c) Id} 

Y, w i cn ' "™ l d n * "" ! (46) 
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The kinetic term acts instead on a state \C) by mapping it to a state \C) (with coef- 
ficient — 1) such that configurations C and C differ by a single plaquette flip. Let us 
now define the set of local rearrangements needed to construct the reduced Hamil- 
tonian fine-tuned to its RK point as the set of plaquette flips (i.e., reversal of the 
arrows on the edges of a plaquette) for all sites i in the square lattice, and all the 
possible initial arrow orientations along the edges. An example of such plaquette 
rearrangements is given in Fig. 2. For reasons that will become clear in a mo- 
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Fig. 2. Example of a decorated and oriented plaquette in the eight- vertex model, and of the 
corresponding local rearrangement generated by the kinetic term in Eq. (44). Notice that 
only the directions of the arrows along the square plaquette (thick line) get reversed. 

ment, let us label these plaquette rearrangements by the set C of all decorated 
and oriented plaquettes t, defined by the site i of their lower-left corner, by the 
orientations of the arrows along the edges of the plaquette, and by the orientations 
of the arrows on all eight dangling bonds that connect the sites of the plaquette to 
(nearest-neighboring) sites not belonging to the plaquette (see Fig. 2). We will also 
denote by £ the (unique) decorated and oriented plaquette that is obtained from i 
by reversing all the four arrows along the edges of the plaquette. Notice that the 
number of elements in C is 16 times larger than the total number of simple ori- 
ented plaquettes (i.e., where only the directions of the arrows along the edges are 
specified). This redundancy, although superfluous for the kinetic operator alone, 
plays a crucial role when writing the potential energy term in the reduced form. We 
can now reinterpret Eq. (44) as 



h smf = i 5>, 



e -Ke e /2p + e +Ke e /2p 



T - + T- 

± Li ~ ± I. 



(47a) 



where 13211 



P t \C) 



and 



T e j\C) 



\C), if£eC, 
0, otherwise, 

\C), ifieC, 

0, otherwise. 



(47b) 



(47c) 
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Here, w e is positive, K is real-valued, Ke t /2 = —{nf 1 — n^)ln(c) — {rig — 
nf^) ln(d) lEiili^l . and C is the (unique) configuration that is obtained from C af- 
ter performing the plaquette flip £ — > £. We stress that the role of the label £ E C 
in Eq. (47a) is two-fold. This label must specify the pair of classical configurations 
that "resonate" through a quantum tunneling process. This label must also carry the 
information needed to construct the difference between the classical configuration 
energies of the two classical configurations that "resonate". For the eight-vortex 
model with no classical interactions between the vertices, i is obtained by decorat- 
ing the oriented plaquettes with the additional information on the orientations of 
the dangling bonds. 

We have reinterpreted the quantum eight- vertex model fine-tuned to its RK point, 
introduced in Ref. W26I1 . using the formalism presented in the past sections through 
Eq. (47a). It is now straightforward to generalize Eq. (47a), e.g., by adding inter- 
actions between the vertices as well as external fields. All this can easily be ac- 
commodated at the classical level, and thus at the quantum level as well, through a 
redefinition of e t , and possibly of C ' . Notice in fact that the definition of decorated 
and oriented plaquettes C r given above, while sufficient to describe an applied 
field or chemical potential on the vertices (as in the original case by Ardonne et 
ah), needs to be modified if interactions between vertices are introduced. For ex- 
ample, in case of nearest-neighbor interactions one has to encode in the definition 
of a decorated and oriented loop I the additional information about the types of 
nearest-neighboring vertices with respect to the ones belonging to the square pla- 
quette, as illustrated in Fig. 3. In the presence of interactions between vertices, the 
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Fig. 3. Example of the change in the definition of a decorated and oriented loop I re- 
quired by the presence of nearest-neighbor interactions between the vertices. The eight 
nearest-neighboring vertices are highlighted with circles. 

coupling constant K can be interpreted as the inverse temperature of the classical 
problem defined by the GS. In addition to changing C by extending the range of 
the decoration of an oriented plaquette, local rearrangements allowing moves in 
classical phase space different than single plaquette flips can be introduced. Any 
combination of these generalizations can change the zero-temperature phase dia- 
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gram of the quantum model. Finally, we stress that the choice for the parameters 
w i can affect in a crucial way the classical stochastics (excitation spectrum) of the 
classical (quantum) system. 

5.3 The quantum three-coloring model 

So far we have illustrated how the formalism of the past sections reproduces quan- 
tum Hamiltonians fine-tuned to their RK points that have already been discussed in 
the literature, and how they can be generalized without spoiling the RK conditions 
(I- III). We devote the following discussion to the construction of a quantum Hamil- 
tonian from an SMF decomposition whose classical counterpart exhibits dynamical 
glassiness. 

Consider the classical three-coloring model on the honeycomb lattice with nearest- 
neighbor interactions ll3vL I38L f39L Eoll . The model consists of classical degrees of 
freedom living on the nearest-neighbor bonds of a honeycomb lattice. Periodic 
boundary conditions are assumed. Any of these classical degrees of freedom can 
assume three different values or colors (say, A, B, and C). The system is also sub- 
jected to the hard constraint that each and every color emanates from any site of the 
honeycomb lattice. The classical configuration space is denoted by S. To each col- 
oring of the honeycomb lattice C E S one can associate an Ising spin configuration 
with the Ising spins defined on the sites of the honeycomb lattice, according to the 
rule that a site is occupied by an up (down) Ising spin if the parity of the sequence 
of colors emanating from this site is even (odd) when encircling the site counter- 
clockwise, say. The classical energy of any coloring C ~ {af\ of the honeycomb 
lattice is taken to be the Heisenberg energy 

E c = -Jj2°> z y (48a) 

Here, denotes an oriented pair of nearest-neighbor sites on the honeycomb 
lattice and J £ M is the spin stiffness. The classical partition function of the three- 
coloring model is then given by 

z(k) = j2 e ~ KEc - ( 48b ) 

The color constraint becomes the requirement that the sum over all the Ising spins 
around any elementary plaquette of the honeycomb lattice, a hexagon, be ±6, in 
the Ising spin representation. We refer the reader to Ref. 114011 and references therein 
for a detailed analysis of the properties on the three-coloring model (48) in and 
out of thermal equilibrium. Suffice here to say that, in addition to the usual Ising 
ordered phases (Neel and ferromagnetic), there exists a large-temperature critical 
phase for < T/J < [T / J) cr in thermal equilibrium 114 lL 14011 . This critical phase 
is expected to terminate in a first-order phase transition to the fully-magnetized 
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ferromagnetic state for sufficiently positive J. The ferromagnetic phase has very 
interesting properties out of thermal equilibrium. At and beyond this first order 
phase transition, the system encounters a dynamical obstruction to equilibration 
(classical dynamical glassiness) that gives origin to a supercooled liquid phase. 
The supercooled liquid phase freezes at a temperature (T/J) < (T / ' J) cr into a 
polycrystallized phase with no interstitial liquid left. In view of the intimate con- 
nection between classical systems out of equilibrium and quantum Hamiltonians 
that admit SMF representations, it is then natural to ask if there is a signature of 
the classical dynamical glassiness in such quantum systems. The detailed answer 
to this question is beyond the scope of this paper. Below, we limit ourselves to 
the construction of an SMF decomposition of a quantum Hamiltonian that has the 
three-coloring model in thermal equilibrium emerging from its GS. 

We define first the separable Hilbert space H of the quantum problem as the span 
of the preferred, orthonormal basis states \C), labeled by the configurations C E S 
of the three-coloring model. Since we are given the classical energy E c of any such 
configuration C, we only need to construct the quantum kinetic energy to complete 
the definition of a quantum Hamiltonian that admits an SMF representation with 
the three-coloring model emerging from its GS. To this end, we note that the only 
classical updates allowed by the constraint involve exchanging the colors along 
closed, non-self-intersecting loops given by the sequence of two alternating colors 
starting at some site i. In the spin language, any such update amounts to nothing 
but flipping all the spins along the loop, i.e., it can be thought of as being induced 
perturbatively by a transverse Ising field H n oc Y,i af . Next, we need to construct 
the set C = {£} of appropriately decorated loops so that: (i) they label all possi- 
ble loop updates, and (ii) they carry enough information for the energy difference 
E c , — E c between two resonating configurations C and C to be written as a function 
s t Both requirements are met when a decorated loop stands for: (1) The backbone 
structure of the loop, i.e., the sequence of sites visited by the loop. (2) The two-color 
sequence covering the loop. (3) The values of the spins that do not belong to the 
loop, but are connected to it by a single (dangling) bond. An example of such dec- 
orated loops is given in Fig. 4. So far, the three-coloring model has been equipped 
with classical stochastics through Monte Carlo simulations using the Metropolis 
algorithm Correspondingly, we choose to implement the classical stochastics 
through a Master equation with transition matrix of the Metropolis form. We are 
now in the position to write the SMF decomposition of the quantum Hamiltonian 
whose GS normalization is nothing but the classical partition function of the clas- 
sical three-coloring model, 

Hsmf =4l>* [e-^Pt + e +K ^ 2 Pj - [f Q + %)] . (49) 

The rescaling coefficient K G M can be chosen at will. The loop energy e g is given 
by the integrability condition (14a): e £ = E c , — E c . The positive coefficient w e 
is chosen to obey the Metropolis condition (28a). The operators P e and T £ j are 
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Fig. 4. Example of a decorated loop in the three-coloring model, as required by the pres- 
ence of nearest-neighbor interactions between the spins, and of the corresponding rear- 
rangement. Notice the information encoded in the loop: (1) the backbone structure, given 
by the sites of all the vertices visited by the AB loop; (2) the two-color sequence covering, 
given by the specific labeling of all the bonds belonging to the loop; and (3) the values of 
the nearest-neighboring spins connected to the loop via the dangling C bonds. 

defined by their action on the basis states \C) 13211 



Pt\C) 



\C), 



and 



\C), 



ill G C, 



otherwise, 



iff G C, 



(50) 



(51) 



0, otherwise. 

As it should be, the GS wavefunction in the preferred basis of the Hilbert space is 

1 



SMF 



Z( K) £ 



-KE c /2 



\c). 



(52) 



It is argued in Ref. that a signature of dynamical glassiness can indeed be 
found in this quantum system when it is brought out of equilibrium through the 
local coupling to a heat bath. 

We shall now refine what is meant above by locality. At the classical level in thermal 
equilibrium defined by Eqs. (48), locality is implemented by a pairwise interaction 
between the Ising spins that is short-ranged. Changing the range of the interaction 
leaves the backbone of a decorated loop £ G C unchanged while it does change its 
decoration. For example, if we modify (48a) by allowing a next-nearest neighbor 
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Heisenberg interaction, we must then decorate a backbone with first and second 
nearest-neighbor dangling bonds. At the classical level out of equilibrium, i.e., at 
the quantum level, locality is implemented by the choice of the dependence on 
i of w e . Since loops of all sizes participate to the classical stochastics (quantum 
dynamics) and since loops involve a number of Ising spins proportional to their 
perimeter, large loop updates (tunneling processes) must be penalized compared 
to small loop updates. Stochastic (dynamical) locality is thus achieved when w £ 
decreases with the loop perimeter in an exponential fashion, say. We conclude that, 
in the Ising spin language, the (weak) coupling to the bath is local if it is between 
the of and a bath of harmonic oscillators. 

We close this discussion of the quantum three-coloring model that admits an SMF 
decomposition by illustrating how it can be thought of as a special point in parame- 
ter space of a more general quantum Hamiltonian, in the same preferred basis. Start 
from the quantum Ising model in a transverse field 



H 



Ising 



hex 



COS [ 27T J2 ^/3 

iGhcx 



(53) 



defined on the honeycomb lattice. The large U limit (i.e., U / J, U/T — > 00) restricts 
the physical Hilbert space to the one obtained by imposing on every elementary 
plaquette (hexagon) of the honeycomb lattice the hard constraint 



cos [ 27T £ cr z /3 

ighcx 



1. 



(54) 



In the same way as the t — J model follows from the Hubbard model in the large 
U limit, the quantum three-coloring Hamiltonian 



= \H [a A + afi - lt (f ere + f- t 



(55) 



follows from (53) when choosing the data (aig,jg) properly II42I1 . The Hamilto- 
nian (49) is another point in the parameter space (a e , 7^) G M 2 specified by condi- 
tions (8) and conditions (14). 

As a second example, we start from the medial lattice of the honeycomb lattice. 
This is the Kagome lattice obtained from the mid-points of the nearest-neighbor 
bonds of the honeycomb lattice. If we attach to any site of the Kagome lattice 
a quantum spin carrying angular momentum one, this gives a three-dimensional 
Hilbert space per site, and a global Hilbert space H of dimension 3^, with iV the 
total number of sites on the Kagome lattice. The three quantum numbers allowed 
for the component of the spin along the quantization (z) axis can be identified with 
the three colors A, B, and C of the quantum three-coloring model. Next, we impose 
two constraints. First, we demand that the coarse-grained magnetization along the 
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quantization axis vanishes. Here, the coarse-grained magnetization is defined lo- 
cally on the honeycomb lattice by attaching to any site of the honeycomb lattice 
the arithmetic average over the spin quantum numbers along the quantization axis 
defined on the bonds meeting at this site. Second, we remove from TL any state 
that has all three spins with vanishing projection along the quantization axis when 
these three spins define a vertex of the honeycomb lattice. The combination of such 
two conditions on H is equivalent to imposing the three-color constraint. Having 
defined the physical Hilbert space H phys , we can then proceed as with the quantum 
Ising spin in a transverse magnetic field and identify a quantum Hamiltonian with 
nearest-neighbor spin-one interaction, say, with Hamiltonian (55). Closely related 
quantum spin-one models can be found in Ref. Il43ll . 



6 Conclusions 

The study of quantum Hamiltonians that are of the Rokhsar-Kivelson type when 
represented in a preferred basis has proven fruitful in the study of quantum critical 
points. This is so because of an intimate relation between the matrix representa- 
tions of these quantum Hamiltonians in preferred bases and some combinatorial 
problems of classical statistical physics. 

We have shown in this paper that the examples of quantum Hamiltonians encoun- 
tered this far, when fine-tuned to their RK points, belong to a broader class of real, 
symmetric, and irreducible matrices that are in one-to-one correspondence with 
matrices known in the mathematic literature as stochastic matrices. Any stochastic 
matrix can be used to define a Master equation of the matrix type that encodes the 
approach to thermal equilibrium of a classical system coupled to a heat bath. It is 
then natural to generalize the known examples of quantum Hamiltonians that are 
represented by matrices obeying the RK conditions (I- III) to any quantum Hamil- 
tonian that admits a matrix representation in this broader class, that we dub the 
Stochastic Matrix Form (SMF) decomposition. Thus, some (not necessarily all) 
quantum phase transitions induced by tuning the coupling constants, say K, [L, B, 
etc., of a quantum Hamiltonian represented by a matrix which is SMF decompos- 
able can be understood as classical phase transitions obtained by tuning the reduced 
inverse temperature K, the chemical potential /i, the external field B, etc., that en- 
ter a classical partition function as intensive thermal variables. This correspondence 
between the quantum GS and a classical partition function extends to the excitations 
of the quantum system and the relaxation modes of the classical system when the 
latter is properly coupled to a thermal bath. 

The correspondence between a quantum Hamiltonian represented by a matrix which 
is SMF decomposable and Master equations of the Matrix type opens the intrigu- 
ing possibility to find a signature in a quantum system for exotic properties of a 
stochastic classical system such as aging and dynamical glassiness. From a more 
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practical point of view, this also suggests a classical Monte Carlo alternative to a 
quantum Monte Carlo simulation when probing numerically the excitation spec- 
trum of a quantum system that allows for an SMF decomposition. 

Finally, we have also proven that any quantum Hamiltonian, up to a shift of its 
eigenvalues, admits a continuous manifold of bases in which its representation is 
SMF decomposable. In this sense, the correspondence between a quantum Hamil- 
tonian construed as an abstract operator acting on some separable Hilbert space and 
the stochastic dynamics describing the relaxation to thermal equilibrium of discrete 
classical statistical systems coupled to heat baths is one-to-continuously many. 

We are indebted to Christopher Henley for a critical reading of our manuscript and 
for challenging us to explore the behavior of the SMF decomposition under a basis 
transformation. 
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